A novel multi-scroll chaotic generator: Analysis, simulation, and implementation
Zhang Gui-Tao, Wang Fa-Qiang
State Key Laboratory of Electrical Insulation and Power Equipment, School of Electrical Engineering, Xi’an Jiaotong University, Xi’an 710049, China

 

† Corresponding author. E-mail: faqwang@mail.xjtu.edu.cn

Abstract

Based on a new three-dimensional autonomous linear system and designing a specific form of saturated function series and a sign function with two variables of system, which are employed to increase saddle-focus equilibrium points with index 2, a novel multi-scroll chaotic system is proposed and its typical dynamical characteristics including bifurcation diagram, Poincare map, and the stability of equilibrium points are analyzed. The hardware circuit is designed and the experimental results are presented for confirmation.

1. Introduction

Lorenz[7] found the first classical chaotic system in 1963, which is a cornerstone laid for scholars to study chaos and its applications. Chaos, which is the complex motion of nonlinear dynamical system, has been widely found and accepted in biology,[1] chemistry,[2] engineering,[3] social sciences,[4] and other fields of science.[5,6] Especially, in the field of electronic circuits, since the typical Chua circuit[8] was presented in 1984, it has quickly become a paradigm to research chaos and a bridge between electronic circuits and the chaos theory.[911] It has a sound engineering application background and lays the foundation for the development of a multi-scroll and multi-wing chaotic system. In the last few decades, multi-scroll or multi-wing chaotic systems have been proposed and some kinds of topologies and shapes of the multi-scroll or multi-wing chaotic attractors have been reported in the literature.[1226] For example, Suykens et al. introduced a quasilinear approach to generating n-double scrolls (n = 1, 2, 3, …).[16] Yalcin et al. presented families of scroll grid attractors including one-dimensional (1D)-grid scroll, two-dimensional (2D)-grid scroll, three-dimensional (3D)-grid scroll attractors with step functions.[21] Tang et al. employed a sine function for generating n-scroll attractors in a simple circuit.[18]et al. presented a switching manifold approach to generating chaotic attractors with multiple merged basins of attraction.[25] Also, Lü et al. introduced the hysteresis and saturated functions series methods to generate multi-scroll chaotic attractors with rigorous mathematical proof and physical realization for the chaotic behaviors.[22,23] Yu et al. presented a general jerk circuit for creating various types of multi-scroll chaotic attractors.[24] As indicated in Ref. [20], compared with conventional simple chaotic system, the multi-scroll or multi-wing chaotic system has a complicated topological structure, so that it is a good candidate for fingerprint images encryption,[27] digital secure communications,[28] entropy of a random number generator,[29] etc.[30,31] Therefore, it is very important and valuable to explore the novel multi-scroll or multi-wing chaotic systems. In this paper, based on a new simple three-dimensional linear system and designing a specific form of saturated function series and a step function with two variables of system, which are employed to increase saddle-focus equilibrium points with index 2, a novel multi-scroll chaotic system is constructed. Both simulation results and hardware circuit implementation are presented for confirmation.

The rest of this paper is organized as follows. In Section 2, a novel multi-scroll chaotic system is presented. The dynamical characteristics of the novel multi-scroll chaotic system are investigated in Section 3. In Section 4, the hardware circuit is designed, and the experimental results are in good agreement with the numerical simulation. Finally, some conclusions are drawn from the present study in Section 5.

2. A novel multi-scroll chaotic system

As is well known, a chaotic system which is a nonlinear system in nature possesses local instability around one or more equilibrium points. This local instability can be identified by linearizing the chaotic system around its equilibrium point to obtain its corresponding linear system, which is easy to judge its stability. Accordingly, this linear system is a key part for constructing the chaotic system.

A novel simple 3D linear system is described by where a is the parameter of system (1). Apparently, is the origin of system (1) and its characteristic equation is given by

According to the Routh–Hurwitz criterion, when , equation (2) has one negative real root and two complex conjugate roots with a positive real part. It is demonstrated that system (2) becomes unstable when a satisfies the condition of and is a saddle-foci point of index 2.

In order to generate multi-scroll chaotic attractors, the main design idea is to introduce additional saddle-focus points by adding some additional breakpoints in the proposed system. Here, a novel multi-scroll chaotic system based on the above simple three-dimensional linear system is constructed by introducing a saturated function series and a special step function with two variables of system.

Accordingly, the novel multi-scroll chaotic system can be expressed as follows: where a, b are parameters and x, y, z are variables of the system (3). In Eq. (4), is one form of saturated function series, K is a positive integer, , , and . Therefore, the saturated slope and turning points of the saturated function series are and (k = 0, ±1, ±2, …, ±K), respectively. For K = 3, the saturated function series of Eq. (4) is shown in Fig. 1. In Eq. (5), can be described as follows:

Fig. 1. Saturated function series when K = 3.

System (3) can generate -scroll chaotic attractors when is selected. Setting the initial condition (x0, y0, , the simulation results from Matlab software can be obtained and shown in Figs. 2 and 3. One can see that system (3) can generate 4-scroll, 8-scroll, and 24-scroll chaotic attractors, when K = 0, K = 1, and K = 5 are selected. Hence, it is effective to use the above method to generate multi-scroll chaotic attractors.

Fig. 2. (color online) 4-scroll chaotic attractors for A = 1.55, α = 0.03, a = 303.03, b = 4.18, K = 0: (a) x and y; (b) xy.
Fig. 3. (color online) -scroll chaotic attractors for A = 1, α = 0.02, b = 1.58: (a) 8-scroll attractor (a = 183.03, K = 1); (b) 24-scroll attractor (a = 190.03, K = 5).
3. Dynamical characteristics of novel multi-scroll chaotic system

In this section, the dynamical characteristics of the novel multi-scroll chaotic system (3) including bifurcation diagram, Poincare map, and equilibrium points are illustrated.

3.1. Bifurcation diagram and Poincare map

Let A = 1.55, α = 0.03, b = 4.18, K = 0 and the initial condition (x0, y0, . The bifurcation diagram which corresponds to the state variable y changing with parameter a is shown in Fig. 4(a). For , one can see that as a increases, bifurcation occurs at the point a = 129.44. In Fig. 4(b), a Poincare map is also given to further observe the dynamical behaviors of the system. The Poincare map is rendered as patches of dense points, and shows that the system is chaotic.

Fig. 4. (color online) Dynamical properties of system (3): (a) bifurcation diagram; (b) Poincare map on the xy phase plane with z = 0, a = 303.03.
3.2. Equilibrium points and its stability

Let the terms on the right-hand side of system (3) be equal to zero, then the following equations will be obtained:

Then, z = 0, , and . Clearly, system (3) has equilibrium points. Here, it should be pointed out that there are two kinds of equilibrium points: and . At the equilibrium points , the Jacobian matrix of system (3) is given by where and the derivative of Eq. (4) at is evaluated by

If , then ; else if , then . Linearization of the system around yields the following characteristic equation:

For , then . Equation (10) can be reduced to For , equation (11) has one negative real root and two complex conjugate roots with positive real part so that all equilibrium points ES are unstable saddle-foci.

For OS , , ±1, ±2, , , then . Equation (10) is similar to Eq. (2). As a result, when , equation (10) has one negative real root and two complex conjugate roots with positive real part so that all equilibrium points OS are also unstable saddle-foci.

Let A = 1.55, α = 0.03, b = 4.18, K = 0, and a = 303.03, then system (3) will have only 4 equilibrium points being OS and three roots of OS will be , –246.26i, and ; let A = 1.55, α = 0.03, b = 4.18, K = 0, and a = −303.03 ( ), then system (3) will have only 4 equilibrium points being ES and three roots of ES will have , , and . In these two cases, although all the equilibrium points are saddle points of index 2, the scrolls of the chaotic attractors can emerge only from OS (when a = 303.03). So it can be deduced that the values of OS are real while the values of ES are virtual in the system. Hence, for , the -scroll chaotic attractors can be generated.

4. Circuit design and experimental results

According to Eq. (3) and Fig. 2, the hardware circuit is designed to have three parts: a main circuit as shown in Fig. 5(c), which has two input ports , and three output ports ( , y, ); a circuit for obtaining function (Fig. 5(a)), which has two input ports , and one output port ; a circuit for obtaining function (Fig. 5(b)), which has one input port (y) and one output port . All the operational amplifiers used in the circuits are all AD711 whose power supplies are ±17 V. It should be noted that (i) when the power supplies of AD711 are ±17 V, its saturation output voltage measured in the experimental setup is ±17 V; (ii) in Fig. 5(b),

Fig. 5. Experimental circuit: (a) circuit for obtaining function ; (b) circuit for obtaining function ; (c) main circuit.

Connect the input and output ports of the designed circuits, then the following equations will be obtained: So Taking , , , , , , , , , , the experimental results can be observed from the digital oscilloscope GDS3254 and shown in Fig. 6. Comparing Fig. 6 with Fig. 2, one can see that they are in good agreement with each other.

Fig. 6. (color online) Experimental results observed from practical circuits for generating 4-scroll chaotic attractor with A = 1.55, α = 0.03, a = 303.03, b = 4.18, and K = 0: (a) time domain diagram of tx and ty; (b) the phase diagram of system (3) on the xy phase plane.
5. Conclusions

In this paper, based on a new 3D autonomous linear system, a novel -scroll chaotic attractor is constructed by introducing a specific form of saturated function series and a step function with two variables of system to increase saddle-focus equilibrium points with index 2. The theoretical analysis, Matlab simulations and circuit experiments are in good agreement with each other and all of them show that this novel multi-scroll chaotic system can generate -scroll chaotic attractors. What is more important is that other novel multi-scroll chaotic topologies can be easily obtained by changing the relationship between the two variables of the system. For example, if which is used here are changed into , or other relationships, one can obtain other quite different novel multi-scroll chaotic attractors.

Reference
[1] Misteli T 2008 Nature 456 333
[2] Cramer J A Booksh K S 2006 J. Chemom. 20 447
[3] Hirata Y Oku M Aihara K 2012 Chaos 22 047511
[4] Murphy P 1996 Public Relat. Rev. 22 95
[5] Kerr R A 2007 Science 317 449
[6] Oxley L George D A R 2007 Environ Model Softw 22 580
[7] Lorenz E 1963 J Atmos Sci 20 130
[8] Matsumoto T 1984 IEEE Trans Circuits Syst. 31 1055
[9] Chen G Ueta T 1999 Int. J. Bifurc. Chaos 9 1465
[10] J Chen G 2002 Int. J. Bifurc. Chaos 12 659
[11] J Chen G Chen D Celikovsky S 2002 Int. J. Bifurc. Chaos 12 2917
[12] Wang C H Xu H Yu F 2013 Int. J. Bifurc. Chaos 23 1350022
[13] Zhang C X Yu S M 2010 Phys. Lett. 374 3029
[14] J H Murali K Sinha S Leung H Aziz-Alaoui M A 2008 Phys. Lett. 372 3234
[15] Wang X Y Lin D Wang Z J 2009 Phys. Lett. 23 2021
[16] Suykens J A K Vandewalle J 1993 IEEE Trans. Circuits Syst. 40 861
[17] Luo X H Tu Z W Liu X R Cai C Liang Y L Gong P 2010 Chin. Phys. 19 070510
[18] Tang W K S Zhong G Q Chen G Man K F 2001 IEEE Trans. Circuits Syst. 48 1369
[19] J Chen G 2006 Int. J. Bifurc. Chaos 16 775
[20] Wang Z H Qi G Y Sun Y X van Wyk B J van Wyk M A 2010 Nonlinear Dyn. 60 443
[21] Yalcin M E Suykens J A K Vandewalle J Ozoguz S 2002 Int. J. Bifurc. Chaos 12 23
[22] J H Chen G R Yu X H Leung H 2004 IEEE Trans. Circuits Syst. 51 2476
[23] J H Yu S M Leung H Chen G R 2006 IEEE Trans. Circuits Syst. 53 149
[24] Yu S M J H Leung H Chen G R 2005 IEEE Trans. Circuits Syst. 52 1459
[25] J H Yu X H Chen G R 2003 IEEE Trans. Circuits Syst. 50 198
[26] Wang F Q Liu C X 2007 Chin. Phys. 16 1009
[27] Han F Hu J Yu X Wang Y 2007 Appl. Math. Comput. 185 931
[28] Tang K S Man K F Chen G 2001 Proc. IEEE Int. Symp. Circ. Syst. 3 787
[29] Yalcin M E 2007 Int. J. Bifurc. Chaos 17 4471
[30] Gamez-Guzman L Cruz-Hernandez C Lopez-Gutierrez R Garcia-Guerrero E E 2009 Commun. Nonlinear. Sci. Numer Simul. 14 2765
[31] Yalcin M E Suykens J A K Vandewalle J 2004 Proc. IEEE Int. Symp. Circ. Syst. 4 581